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Abstract 

A dual Ginzburg-Landau model of the knotted chromo-electric flux-tube is revisited, 
in which the covariant decomposition of gluon field and the random phase approximation 
are used, ft is shown that the SU(2) QCD vacuum is of type-II superconductor, with the 
Ginzburg-Landau parameter k = \fi, being independent of the magnetic condensate and 
strong coupling used, and consistent with the lattice data. The mass spectrum of a number 
of /j meson states with J < 2, which are taken to be of glue dominate, are computed 
with help of the energies of the knotted (linked) QCD fluxtubes. The low-lying /j states 
[leqX.lGeV) are shown to be associated with the string excitations of the types 2^ and 
X = 3]^ in knot topologies. 

PACS number(s): 12.38.-t, 11.15.Tk, 12.38.Aw 

Key Words QCD vacuum, Ginzburg-Landau model, Glueballs, / meson states, knotted 
fluxtube 



1 Introduction 

According to the theory of strong interaction, quantum chromodynamics(QCD), there will be 
pure-glue excitation known as glueball, and the qq states with explicit glue (known as hybrid 
mesons), which can have quantum numbers forbidden to the qq systems in the naive quark 
model. These gluonic excitations have been the subject of the extensive experimental and 
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theoretical studies(see, for instance, [H [2]) as the identification and observation of these objects 
can be a good test of QCD at the low-energy limit. After the early work of the bag-modcl[3j 
of the glueballs, with mass prediction about 1.0-2.0GeV/c^ , many glueball (or, glue-rich meson) 
states were explored with various approaches, including QCD on the lattice, the QCD sum rule, 
the constituent gluon models, and etc., leading to a deep relation between the properties of these 
states and the structure of the QCD vacuum[Tl H]. Owing to the mixing with the normal qq 
states, to accommodate these gluonic excitations (except for the lightest tensor sector) remains 
to be an open question[ll El |3]. 

The present work revisits the low- lying fj meson states (listed in the Particle Data Group [5]) 
using the dual Ginzburg-Landau(GL) model proposed in our previous work[6], which is based 
on the reformulated Yang-Mill(YM) theory in terms of the field decomposition[71 El E] of gluon 
variables. We show that the vacuum of the SU(2) QCD is of a type-II superconductor type 
and the GL parameter is k = \/3 for the vacuum condensate, being consistent with the lattice 
simulation [TUj and independent of the magnetic condensate as well as the strong coupling that 
are used in the calculation. We further compute the mass spectrum for a number of the 
low-lying fj meson states assuming that the these mesons are mainly formed by the tightly 
knotted(or linked) QCD fluxtubes, that is, the knotted(or linked) QCD strings (Hereafter, we 
will also use the term knots for the closed types of strings including links). The random vacuum 
approximation is used for the vacuum and main features of a few low- lying fj states, the 
candidates for the glue-dominated states, are examined in associated with the knot topology of 
the chromo-electric fluxtubes for J < 2. Though the starting point for model construction is 
the two-color gluodynamics, the relevance of the present work to the real QCD can be inferred 
from the fact that A^-color QCD can be expressed as a sum over copies of two-color QCD [11] as 
long as the Abelian dominance is valid for QCD[T2], as shown in lattice QCD[TH 1 114 t I15j. 

The idea of the gluonic fluxtubes [HI [T7] can be traced back to the early days of QCD when 
the resonance string(i.e., the QCD string) was used to describe the confining force connecting 
the quarks [18]. The formation of a gluonic fluxtube between two widely separated quarks is 
widely accepted, and is supported by the lattice QCD simulations [131 dU [T5]. Meanwhile, it is 
expected that the glueballs, if exist, are intimately related with the closed fluxtubes that can be 
omitted, for instance, by a long linear string between the quarks. Using a flux-tube model, Isgur 
and Paton[19] predicted that the lightest glueball has the quantum number J^*^ = 0"'"'" and a 
mass of about 1.520Gey/c^. The further predictions made later by the flux-tube model[20] for 
the three lightest glueball masses are consistent with the lattice calculations [Tl \n] . Assuming 
the closed (knotted) configurations for the gluonic fiuxtubes, the glueballs were explored by 
using the various actions, such as the Nambu-Goto action in the case of the circular string [22], 
the nonlinear sigma actions [23] and its extensions |24j in the case of knotted strings. For the 
early attempts to model hadrons in terms of closed strings of quantized fluxtube, see Ref . [25] . 
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Our study is motivated by the ideal representation [26] of the knot geometry, which provides 
a model-independent relationship between the length-to-diameter ratio and the average crossing 
number of knot with the given topology. This representation complements the nonlinear dy- 
namics of the stable knot found in field theory [271 EH], ^-nd is greatly helpful in evaluating the 
glueball-like meson spectrum in terms of the knotted objects in a bag model of fluxtube[29j. The 
calculations in this work involve the Nielsen-Olesen(NO) vortices in the dual GL theory as an 
effective description of the gluonic fluxtubes, for which the knot geometry previously explored 
was utilized to calculate the energies of the low-lying glueball-like meson states, similar to that 
in [29]. 

2 The order parameters in long-distance gluodynamics 

We begin with the reformulations of SU (2) gluodynamics via the covariant field decomposition of 
SU(2) gluon variables, known as Cho-Faddeev-Niemi(CFN) decomposition [3 [3 IS]- The gluon 
field (the arrow denotes the three color indices a = 1,2,3, along the generators r") is 
decomposed into [7] = A^n + C^ + X^,, in which = A^-n is an Abelian gluonic potential, 
n{x) an unit isotriplet in color space, (7^ := g~^d^n x h the (non-Abelian) magnetic potential, 
and Xfj^ (normal to n) an covariant field, having to be constrained by two extra conditions [30j. 
We use the simplest choice of the condition[9j, = g'^^cpid^h + g~'^4'2d^n x n. With this 
change of variables, the YM Lagrangian becomes 

^YM ^ _ 1 [^^^ _ ^H^^f - ^{{n^u - ^//;..)(V^V)V> + /i.c}, (1) 

where F^,^ = df^A^ - d^A^, Z{(j)) = 1 - <f) = (fyi + i(j)2, n^i, = rifj,y{dnf - d^h ■ dyh, 
V^(j) = (<9/i ~ W-^ij)4> is the U{1) covariant derivative induced by the gauge U{1) rotation 
U{an) = exp{ian"'T"') around direction h. The chromo-magnetic field H^u = h ■ {d^h x duh) is 
partially dual to the chromo-electric field Ffj^iy, as shown in |8l[9]. Such a change of variables in 
the reformulated dynamics ([T]) implements the Abelian projection [12] in a covariant way[2]. By 
choosing h{x) in ([1]) as an infrared order parameter and using renormalization group analysis, 
the Skyrme-Faddeev(SF) model[27], which supports the knotted solitons with nonzero Hopf 
charges, was proposed to model the low-energy dynamics of quantum YM theory [9]. 

To see the implications of order parameter fi{x), it is helpful to re-examine the connection of 
the CFN field decomposition with the Abelian projection from the viewpoint of the basis change. 
As shown in [3U] the choice of the constrains on implies to choose a local gauge transformation 
which maps the generic gluon field into a gauge-fixed surface in the space of the gluon fields, 
which yields that the map will necessarily be singular somewhere in spacetime. The magnetic 
degree enters through the topological variable n{x), which provides the knot dynamics [271 ES] 
and classifies the physical region V considered (as a mapping from V to S^) according to the 
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homotopy TT2iV). We note that when the U{1) symmetry (rotation around h by angle a) is 
left unbroken, h serves as the mapping transformation from the asymptotically free gluon Afj_, 
which is represented in terms of the global basis to the infrared variables, which are 

represented in terms of the local basis {n, c?^n,5^rexn}, with nontrivial metric in the gauge group 
space. The QCD vacuum differs from the perturbative one owing to the nontrivial homotopic 
class of the map n(x), or, equivalently, to the singularities (magnetic charges) in the magnetic 
potential C^, namely, the zeroes of the map h{x). The field decomposition in terms of the local 
basis collapses when h becomes globally fixed n{x) — )■ no(i-e., the norm of d^fi vanishes) so that 
{n, d^n,dfj_h x h} degenerates. In the latter situation, one has instead to go back, discontinuously 
in the field mapping, to the usual asymptotically free gluon variables = A°'{x)t"' in the usual 
matrix representation of A^. This discontinuous transition in the local gauge- fixing differs 
the asymptotically free phase of QCD from the confining one. Thus, the nonvanishing vacuum 
expectation value(VEV.) {{dn)"^) is required for the CFN field decomposition to be a true change 
of variables. 

Another order parameter arises by looking the confining phase of the theory ([T]). Taking the 
pure Abelian gauge = so that classically F^^i, = 0, as should be in the presumed vacuum 
condensate, the theory ([1]) becomes 

in which (p = p(x)e**^^^ has been used. One sees here that Z{p = \(j)\) resembles the dia-electric 
factor in the dia-electric soliton model [31] and the gauge-invariant kernel in the effective model 
of confinement j32l if one interprets the background media associated with (p as the QCD vacuum. 
Explicitly, Z{<j) — )• 0) = 1 corresponds to the normal vacuum and Z{p — t- u) 7^ (here, v = {p) 
is positive constant) to the condensate. Indeed, taking the limit cj) ve^^°, the Eq. can be 
further reduced to 

+V{n-h,0 + -- - . 

Here, the second term quadric in derivatives of n and potential terms V{n -h,^) are added based 
on the renormalization group analysis, as done in [9l[23]. As we can see, ([3]) is an extended version 
of the SF model p7] (see also ^^). Given the nontrivial configuration n{x), the knot is classified 
by non-zero Hopf charge Q = l/(327r2) f d^xe^^^CiHjk, with Cj (having no local form in terms 
of n) defined mathematically by dC = H (The notation H = Hj^dx^ A dx^/2 = (n,dn A dh)). 
Assuming there is a localized excitation G of 77,(2;) in subregion Vq C V and using the virial 
theorem, the energy density estimated by ([3]) is about Hq oc 2((9^n)2) -|- Vmim which is in 
consistent with {{d^fif') ^ in Vg, while {{d^h)"^) = far away from Vq due to the finite energy 
condition. The situation is opposite for (j), which tends to zero at the core of Vq and to the 
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nonzero condensate v = (p) outside Vq, as it should be according to the resemblance between 
Z{p = in ^ and the dia-electric factor. 

The quanlitative behavior of fi and (p in G can be related to the chromo-magnetic symmetry 
breaking, where the chromo-magnetic symmetry H^j is defined by [8] 

Df,h{x) = {df, + gAf, X )n(x) = 0. (4) 

Observed that the CFN decomposition implies D^h(x) = giX/^i x n), one sees that the breaking 
of Hm is amount to the nonvanishing VEV. of [D^nf'^ namely, 

/ g\{X^yihf) 

= {{(pidf^h xh- (j)2df,nf) 

= {{4 + ^l){d,hf), (5) 

in which the reparameterization [9j for is used. In the case of weak correlation between the 
off-diagonal variables (0, n), implies 

/ (|<AP), Outside core(yG), .g. 
0^{{d^h)^),inVG, 

which agrees with the existence of the order parameters n and (j). In short, the arising of the 
order parameter h and (p are associated with the breaking of Hm defined by Observed that 
{X^ X n)^ = X'^ in which Xf^ • n = 0, the chromo-magnetic symmetry breaking ([5]) is amount to 
the X-field condensation 

0^{{X,f) = m\g-'d,nf). 

Such type of condensation, referred as the off-diagonal gluon condensation, was seen in the 
lattice simulations [33] for gluodynamics. 



3 The random phase approximation 

To have the effective dynamics of the Abelian-Higgs multiplets ( , 0) , we need the further 
approximation associated with the vacuum of QCD. We will use the random phase approxima- 
tion(RPA), which can be due to the residual U{1) symmetry shared by the dual dynamics ([T]). 
To see this, it is convenient to use two normalized basis {ei, 62}, such that ei x 62 = n, to describe 
the 2-dimensional plane '^(x) which is normal to h{x) at x. One can use a transformation from 
{61,62} to {e+,e_} explicitly for the isomorphic map: SO{2) — U{1). 

Take to be the direction vector of O^n, one can then write, as d^h is in the plane 

6^ = 0^61+02 62, 

= 0/^6+ + a^e_, 
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where are the two real four-vectors such that (a^)^ + (a^)^ = 1, and 



e+ = i(ei + ^62), e_ = i(ei - ^62). 



It is easy to show that 



e+ • e_ = 1/2, e+ • e+ = e_ • e_ = 0, 

e+ X e_ = — = — e_ x e+, 
e+ X n = ie+ , e_ x h = — ze_ , 

The last equation of ([5]) implies that a^{x) = e*^*"*^^), with ^^(x) a real phase. 
From the equations ([8]), one has 

= cos{0^ - By), 

= sm{e^ - Oy), 
e^-ey + i{e^ x e^) ■ n = e'^^^-^-'\ 
Under the small gauge rotation U{ah), with a a small angle, one can show that 

fi/i + (5e^, 
= + X (an), 
= (a^ + aa2)ei + (a^ - aai)e2, 
= a'^e+ + a^e_, 

where 



(8) 



This means that the gauge rotation U (an) corresponds to 

e/^"^^("+^'')e+ + e-("+^'')e_, (9) 

or equivalently, to the phase shift 0^ ^ 9^ + a. 

By writing d^h = M{x)e^ and using ([7]) and ([8]), one has for the SU{2) gluon field 

M(x) 

A^, = A^,h + —-^ [{i + 0)a^e+ + h.c] (10) 
one can show that under U{ah) 

(^^ + ^) n + ^ [(i + + ia0)a^e+ + h.c] , (11) 
'A, + ^)h+MM [(i + e-<A)a^e+ + /i.c] . 
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Comparison (jlOp with (jlip shows that the rotation U{ah) leads to transformation in the variables 
{A^, (j)) as exactly as that for the Abelian Higgs multiplets: — ^4^ + d^a/g, cp e*°(/). This 
indicates that the unbroken U{1) symmetry shared by the dual dynamics ([1]) is the local gauge 
rotation:[/(an). 

We introduce the RPA such that (e^) = 0, which, by noticing the residual symmetry ([9]), 
means 

{e'^^) = 0. (12) 
Moreover, in the case that l^. ^ f one has, from (|12p 

(ei(eM-e.)) ^ (e*^'')(e-^^'^) = 0. (13) 

The RPA introduced here assumes 6^ to be random distributed in the QCD vacuum. As will 
be explored in the following section, the RPA is very useful to find the effective U{1) dynamics 
of the collective variables (yl^,</>). 

In order to describe the knot-like excitations, the SF-like dynamics of the order parameter 
n, similar to ([3]), are usually used, as done in [341 124j . Owing to the difficulty for extracting 
the parameters in (l3|) , an alternative approach is to use the dynamics of (^^ , (/>) to describe the 
fine profile of the chromo-electric fluxtube, and to extract the energy of the closed fluxtube by 
utilizing the values of the universal invariants for knot geometry [26j . which is main propose of 
this paper. 



4 The dual Ginzburg-Landau model 

From the section 2, we know that the order parameter (t){x) in ([1]) is suitable to play the role 
of soliton field interpolating in between the two vacua: 4'{x) = and </>(x) = 0). Writing 
0(x) = <I>(x) + 5(j), where <I>(x) is the complex condensate and 54> its quantum fiuctuation, one 
has a nonzero correlation 

{(t){x)^\y)) ^ $(x)$*(2/),for x° > . (14) 

To find an effective model for [A^,(j)) starting from ([1]), let us take the mean field approxi- 
mation 

df,n{x) = Me^{x), (15) 

with the unit vectors defined in d?]), and M = ((9^n)^)^/^ (no summing over /u) a constant 
in spacetime. The equation (jlSp yields 

{d,hf = M2{(eo)2 - E-=i(e,f } = 

H^, = M^V, (16) 
h^y = n- (e^ X eu) = siwO^y 
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with 6^1, = 0^ — 6v the angle between the directions and in the internal space. In addition, 
one has 

dfj,n ■ dyh = cos ^^jy, 
n^^ = -M2(2r/^^ - cos 6^^), (17) 
H^.H'^'^ = M%,h>^^ = ^ Z^^^il - cos2V)- 

Defining the magnetic condensate H by = {H^,yH'^'^) and using the RPA, one finds, by 
using (HZD and (USD, 

= QM"^, (18) 

namely, = H/^/6. In deriving (fT8|) . we have used J^fiu 1 = 12, and (cos 26*^1,) ~ according 
to the RPA. We see from ()18p that the very existence of the nonvanishing VEV. (((9n)^) implies 
the chromo-magnetic condensation, {H^j^yH^^^) oc M^, which is due to the breaking ([5]) of the 
magnetic symmetry Hm- 

It is important to note here that the RPA ignores the possible nontrivial bending (or, 
twisting) of the h{x) orientation, which is crucial to identify the knot-like excitations. In the 
case that h{x) = hci{x) + 5h, where the variation of the classic part nci{x) is small(namely, 
the contribution to M arises mainly from the quantum fluctuation 5n), one can take hd to be 
approximately parallel. Then, upon using [T^ and the mean field approximation (|15|) as well as 
(fT6|) . the gluodynamics ([1]) becomes 

/Vf2 

+ 4^{[2r?^. + e*^-](V'^0)tV> + h.c}, 
which, after utilizing becomes 

^ -\f',. + -^((V^0)tv» - ^{Z{<Pf). (19) 
Here, the following relations, which are due to the RPA, are used, 

( V) - 0, (/i2 J ~ 6, 
(e»^/--(V0)tV</>) ~ 0. 

Using the Wick theorem and the Bose symmetry of the scalar field, one can show 

((0t</,)2) = + + 

= 2(0t0)2, 

and 

(Z(</>)2) = (1 + (0t^)2 _ 20t^) 

~ 1 + 2(^>*^>)2 - 2$*4> 
= 2(|«>|2- 1/2)2 + 1/4. 
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where ([H]) is applied so that ((V^(?!))'^V^(/)) = (V^<I>(x))* V^<I>(x). Rescahng the scalar <i> to that 
with dimension of mass, 



^9 

we obtain, from ()19p . the effective dual GL model given by 



IDGL ^ _1^2^ + 1(5^ _ ^^^^)^|2 _ _ (20) 



In terms of the re-scaled complex the potential in (|2U|) is 

A; 
4 



y(<I>) = ^(|ci)p-^2)2^ (21) 



with the parameters given by 

A = Vl2g, 

One sees that the effective model ()20p for the gluodynamics takes the form of that for the 
dual superconductor [53]. It is remarkable that the scalar potential (|2ip assumes exactly the 
Mexico-hat form and ensures the dual Meissner effect for confining the chromo-electric field A^, 
since both of ~ and g are positive. We also note that the ensuing dual superconductor 
picture with two vacua <I> = and <I> = ?; in ()20p agrees with the vacuum picture discussed in 
section 2. 

It is known that the dual GL model (|2Up admits the NO vortex solution |16] , with two length 
scales: the coherent length ^ = l/m$ and the penetrating length \l = l/ruA- The mass scales 
m$ for the Higgs-like field <I> and for the chromo-electric field are fixed by the explicit 
form of the potential (|2ip . Writing in terms of A and v in (|22p . or, equivalently of the magnetic 
condensate H, they are 

if, 



ruA = V2gv = ^, 

which depend merely upon H = \/6((5n)^) and are nonzero when the magnetic symmetry 
broken. Given the scales (|23p . one readily finds the GL parameter (defined by k = for 
the condensate vacuum to be 

^ = — = (type-II), (24) 

niA 

which is independent of the magnetic condensate H. The GL parameter given by (|24p predicts 
the vacuum of the SU(2) gluodynamics to be of the type of type-II superconductor, in nicely 
consistent with the lattice data k = 1.702(= 0.16/?7i/0.094/m) for SU(2) gluodynamics [TO]. 
The very fact that the potential (|2ip in (|2Up assumes the Mexico-hat form gives an independent 
argument for supporting the dual superconductor mechanism of the low-energy phase of QCD 
proposed by Nambu, 't Hooft and others [T71IT2]. The type-II superconductor was also confirmed 
in some of the lattice simulations, in which k = 1.04(= 1.3614Gey/1.3123Gey) ^37] and k = 
1.49(= 0.164/m/0.11/m)[36j were predicted. 
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5 Glueball-like mesons as knotted fluxtubes 



Due to the U{1) gauge symmetry for which it2{U{1) = 1), the dual GL model (|20p does not 
allow the stable solution of closed fluxtube being against self-shrinking. To describe the knotted 
gluonic excitations, it entails to have a nonlinear dynamics of knotted configuration, such as SF- 
like model ([3]), which is difficult to solve for now. A way out is to use the Nambu-Goto action [22] 
to prevent the self-shrinking instability of the closed fluxtubes in the GL model. In this section, 
we use (PU|) as a model of the gluonic flux-tube profile and calculate the glueball spectrum by 
taking into account both the string tension and the twisting of the knotted fluxtubes, given that 
the geometric ratios of the length to the diameter are known for a set of knot types [26]. 

Consider a slice of the vortex cross-section, with the cylindrical coordinates (r, 9) in it. We 
search for, as usual, the static NO(or, the Abrikosov) vortex solution to the model ()20p . in the 
Coulomb gauge (V • A = 0), in the form 



gr I 



A{x) = 9A{r) = -'^[l-F{r)],Ao = 0, 



with the boundary condition F(r — )• oo) = 0, p{r — )• oo) = 1. The chromo-electric field is given 
by {E, B), with E = in the present case, and 

B = V xA = z-(^] . (25) 
gr \dr J 



The static energy for (j20|) is 

^DGL^ |d3^|:|! + |(9^_i5A0a>|2 + y($) + g!| (26) 

with V{^) given by ()2ip . In terms of the vortex profiles {F{r),p{r)), the energy (|26p becomes 



E= d^x 



' ^ ' ■ +^ hr +^F^p' + V{p) 



^ \gf J \dr J \dr 

where the constant energy (e oc J H'^d^x) in ()26p is treated as the zero-point energy. Using ()2ip 
and ([2^ . one has for the potential. 

Introducing a dimensionless variable x = r/a, with a a length scale, the static equations of 
motion for (|27p are 

d'^F IdF , 
ax^ X dx x'^ 
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with the controlhng parameters given by 



A 
B -- 



va) 



{Xva/V2f. 

The asymptotic behavior of the vortex profiles F{x) and p{x) reads 

p{x ^ 0) = c'^KniVBx) 



(29) 



and 



F{x ^ oo) = CnxKx{\fAx) ^ CnVxe-^"", 



p{x 



oo 



1 - CLKn(V2Bx) 



1 



CL 



riBx 



(30) 



(31) 



with Cn,c'^,Cn,C^ the constants for given n, and Kj{x) the Bessel function of second type. The 
expressions (|30|) and (|31|) are useful and will be incorporated in the specification of the boundary 
condition in the process of the numerical relaxation for solving ()28p in finite (but large) interval. 

Assuming that the vortex has a finite length L along z direction and integrating the angular (^) 
part, one has, from ([2j 



2-kL 



{n/gf 



jglue ^ ^2jB ^ 



-I 



N 



(32) 



with three dimensionless integrals given by 

jglue _ j-co dx (dF\'^ 
Jo 2a; V da: / ' 



B 



dx 



X 



(33) 



N 



The first term in (|32|) is the chromo-electric energy and roughly scales as ~ L(tr$|)/(7ra2), due 
to the flux conservation, with the flux of chromo-electric field ()25p in the area ira^ of the 
vortex cross-section. Since the three integrals in (j33p are all positive for the NO vortex, the 
stable vortex exists, for which the first and third terms in ()32p balance. 
The minimization of ([32]) with respect to a yields 

^T7 



o(= 



H 



which depends explicitly on the vortex quantum number n. In deriving ([2 
the other hand, using ([22]) . (P^j) and (fM|) . one has 



A 



(34) 

22]) are used. On 

(35) 



This means that the energy depends on a„ not only through (|32p explicitly, but also implicitly 
through the integrals (|33p . which rely on a„ through the controlling parameters A and i3(both 
of them ~ a^). 
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For given input (A,v), or equivalently, (g,H), we numerically fix the the length scale a in 
the case n = 1 ~ 7, in two steps: (1) to find the vortex profile {F, p) by solving the GL equation 
()28p using relaxation method with A and B given by (|29p for a given a; (2) to find the scale 
a = an for which the string tension a = E^/L is minimized by calculating E^/L as a function 
of a through (f32]l and (f33]l for the profile given in (1). The optimal cl and the corresponding k^i^ 
including three integrals (/^'"^^ I^), are listed in Table I. The numerical results are shown in 
FIG.l for the profile {F, p) and FIG.2 for the tension o"(a) as a function of a, for n up to 4. The 
relation ([35]) are found to be fulfilled for calculated a„ in (pl|) . In the process of the numerical 
relaxation, the asymptotic behavior (j30p and ()3ip are incorporated in the boundary condition 
and the initial vortex profiles Finitiai{x) = sec h{^/Ax) and Pinitiai{x) = tanh(\/2-Bx) are used. 





1 


2 


3 


4 


5 


6 


7 




1.0019 


1.2982 


1.5204 


1.5945 


1.7426 


1.9648 


2.2611 


rn{Gev~^) 


0.5012 


1.0490 


1.4235 


1.6156 


1.8551 


2.1420 


2.5520 


jglue 


0.1941 


0.1941 


0.1942 


0.1942 


0.1936 


0.1933 


0.1944 




0.6770 


1.5893 


2.8129 


4.3547 


5.7913 


6.5896 


7.3093 


jN 


0.0370 


0.0418 


0.0461 


0.0490 


0.0525 


0.0601 


0.0687 




1.2850 


1.6650 


1.9500 


2.0450 


2.2350 


2.5200 


2.9000 



Table I 



1 1 1 1 1 1 1 1 1 1 1 




B for n=l 


\ / ' 


Fforn=l 




p fpj- 




F for n=2 




p for n=2 




F for n=3 




p for n=3 




F for n=4 




p j-Qj. jj=4 







0.0 0.2 0.4 0.6 0.8 1.0 

r(fm) 

Figure 1: The vortex profiles for n = 1 ~ 4. F decreases with r and changes slightly for different 
n, while p increases faster as n increases. The plot is also given for the chromo-electric field 
i?(with unit of GeV"^), which decreases from 0.1866^6^^ (at r = 0) to 0(at r = oo). Here, the 
strong coupling and the magnetic condensation are taken to be = 3.2198, H = 0.5246Gey^. 

It can be seen that ()32p describes the linear rising (~ crL) of the confining potential for large 
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0.6 

r(fm) 



Figure 2: The string(vortex) tension a = E'^ / L as a function of the length scale a for n = 1 ~ 4. 
Here, g = 3.2198, H = 0.5246Gey2. 

L. When L decreases, there will be corrections to (j32p due to short-range gluonic interactions. In 
QCD we can expect, quite generally, the full energy for the low- lying spectrum of the confining 
fluxtube with the fixed ends at distance L to be|381 [39] 

E{L) = aL + e+'j- + '^ + 0{^). (36) 

Here, a is the string tension, e the the zero-point energy in (|26p . rmxjL (m = 1,2,---) the 
vibrationally-excited energies of the string flux, and cq = — -7r/12 is the Casimir energy of zero- 
point fluctuations of the string [58]. Combining (p6]) with ([32]). one obtains for the energy of the 
knot-like gluonic fluxtube. 

^2 



En,m = 4:TTe{K)rn 



an 9 



2 



2e(A)r„a„ 



where r„ = i?/a„ is the ratio of the core radius R of the vortex to an in ([M]) . and e{K) = L/ (2i?) 
a topological invariant which is universal for a given knot types [26]. The factor r„ depends merely 
upon n and can be fixed through the numerical relaxation for (j28p . In principle, e{K) can be 
evaluated by the dynamics of knot like ^ . This is not available by now since we fail to fix the 
explicit form of the potential V{n • /i, ^) in ([3]). Fortunately, e{K) were previously determined 
for a set of knot types (K's) via the Monte-Carlo simulations [26], partial of which are listed in 
Table II and HI for our references. 

We use the fluxtube energy (j37|) for these knot types(denoted hy K = 2f , 3i, etc. in topology) 
to model the spectrum of the glueball-like mesons. The computed spectra are shown and 
compared with the data of the partial fj states (J < 2) (also a few rj states) in the Table 
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II in the case of n = 1 and of some low m's. In the evaluation of the glueball-like meson 
energies using (|37|) . the strong coupling and chromo- magnetic condensate are chosen as 

g = 3.2198, H = 0.524QGeV'^ (38) 
so that Aas/S = g'^/{M) = 1.1 (see Ref.pjj) and 

\,yac\=H^/{^g^) = {Q.2AGevf- 
V = ^/Hj2/{^/^g) = O.lOieGeF 

with evac the vacuum energy density, given in the lattice computation |4Uj . H is fixed based 
on ()23p through the lattice data for the coherent length ^ = l/m$ and penetrating length 
Xl = ^/fnA- The numeric result for is oi = 1.2850 for n = 1 with the tension 

a = 0.125Gey^ 

and the other cases n > 2 are excluded since they yield anomalously high string tension, such 
as fj = 0.339, 0.695, etc., see FIG. 2. 

Different with ([38|) . we alternatively choose, by taking v = 0.1094Gey and evac = — O.TGeT^, 
the input values of the parameters which change H slightly, and thereby the scalar condensate 
V = ($) (see ^) so that 

g = 3.2198, H = 0.6080Gey. (40) 

The numerical data corresponding to (j40|) are shown in Table II and III for the mass spectrum 
of a number of low-lying fj states with J < 2. 

It can be seen from Table II and Table III that a remarkable agreement was reached between 
the experimental and predicted spectrum for the meson fj states, showing that the most of 
the meson fj states can be identified as the knotted string excitations of the types (n, m) = 
(l,j),j = 0,1,2. A little different suggestion for these identification are given in Table IV, 
with the prediction for glueballs with quantum number J^^ = , , 0^ , compared to the 
lattice data. 

We account for the results in Table II and III as follows: a number of the meson fj states, 
when taken to be the glueball-like (i.e., the glueball dominate) states, can be viewed as a knot 
excitation of the chromo-electric fiuxtubes with knot types n'^(n = 1 ~ 5, / = 1, 2) in the topology 
terminology,in which the low-lying (< 1. 71 GeT^/c^) glueball-like meson states are best described 
by the fluxtube of the knot types 2\ and 3^, with (n, m) = (1, j), j = 0, 1, 2. It is unknown in our 
framework why some ofvibrational modes lack the experimental counterparts, but it is known 
thatthe knot types K = 2f,3^ corresponds to the low- lying fj states because they have the 
simplest topologies among the knot types 2^,3^, in the sense that they are non-shrinkable and 
they are shortest in the unit of the string diameter. While the minor variation, roughly 2{)QMeV 
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per unit of m, in the spectrum arises from the vibrational modes of the fiuxtube, the spectrum 
of glueball-hke mesons are mainly due to the knot geometry of fiuxtubes. 



States 


Mass" 


m 




eiKY 


En,m{G) 


/o(600) 


400-1200 


1 


2i 


6.2832 


674.6 


/o(980) 


980±10 





2? 


12.6 


976.0 


/ii(1170) 




2 


2i 


6.2832 


1246.9 


ai(1260) 




1 


2f 


12.6 


1262.2 


/o(1370) 


1200-1500 


1 


2? 


12.6 


1262.2 


/o(1500) 


1507±5 





3i 


16.4 


1500.1 


/o(1710) 


1718±6 


1 


3i 


16.4 


1719.4 


ry(1760) 


1760±11 


3 


2? 


12.5664 


1757.0 


r/(2225) 


2220±18 


1 


2f *0i 


20.8496 


2278.4 


/2(1910) 


1915±7 


2 


3i 


16.4 


1938.7 


/2(2010) 


2011+^° 


4 


2? 


12.5664 


2015.5 


/o(2100) 


2103±7 


3 


3i 


16.4000 


2100.4 


/2(2150) 


2156±11 





4i 


21.2000 


2153.6 


/o(2200) 


2189±13 





4f 


21.4000 


2180.8 


/2(2300) 


2297±28 


4 


3i 


16.4000 


2298.5 


/2(2340) 


2339±60 


1 


4? 


21.4000 


2348.8 



Table II 
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States 


Mass" 


m 


K" 


/ T7-\f' 

e{Kf 


T— 1 / 


/o(600) 


400-1200 


1 


2i 


6.2832 


674.6 


/o(980) 


980±10 





2{ 


12.6 


976.0 


/ii(1170) 




2 


2i 


6.2832 


1246.9 


ai(1260) 




1 


2{ 


12.6 


1262.2 


/o(1370) 


1200-1500 


1 


2! 


12.6 


1262.2 


/o(1500) 


1507±5 





3i 


16.4 


1500.1 


/o(1710) 


1718±6 


1 


3i 


16.4 


1719.4 


??(1760) 


1760±11 


3 


2i 


12.5664 


1757.0 


??(2225) 


2220±18 


1 


2| * Oi 


20.8496 


2278.4 


/2(1910) 


1915ib7 


2 


3i 


16.4 


1938.7 


/2(2010) 


2011+6° 


4 


2f 


12.5664 


2015.5 


/o(2100) 


2103±7 


3 


3i 


16.4000 


2100.4 


/2(2150) 


2156±11 





4i 


21.2000 


2153.6 


/o(2200) 


2189±13 





4? 


21.4000 


2180.8 


/2(2300) 


2297±28 


4 


3i 


16.4000 


2298.5 


/2(2340) 


2339±60 


1 


4? 


21.4000 


2348.8 



Table III 



The data from the PDG summary tables and[2]. ^ The knot types in terms of notation n^, 
means a link of I components with n crossing, and occurring in the standard table of links() on 
the /cth place. The data from the simulations [26| except for 2^ and 2f * Oi. The data from 
Lattice simulations. 



States 


Mass'^ 


m 






En,m{G) 


0++ 


1710±130 


1 


3i 


16.4 


1678.5 


0-+ 


2560±155 





5i 


24.2000 


2561.1 


0+- 


4780±290 





92 


40 


4703.3 



Table IV 



6 Summary and discussions 

The dual dynamics of the SU(2) QCD is revisited using the covariant (Cho-Faddeev-Niemi) 
decomposition of the gluon field. Assuming that the chromo-magnetic symmetry is broken and 
using the random phase approximation, we show by deriving a dual Ginzburg-Landau model for 
SU(2) gluodynamics that the QCD vacuum is of type-II superconductor. The Ginzburg-Landau 
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parameter is shown to be k = \/3, which is independent of the chromo-magnetic condensate H 
and the strong couphng Qs used, and agrees well with the lattice simulation. The mass spectrum 
of a number of low-lying fj states are calculated for the chromo-flux quantum n = 1 and found 
to be in a good agreement with the recent lattice data, by adding the energy correction arising 
from the twisting of the knotted (linked) QCD fluxtubes and the zero-point energy to the energy 
of the Nielsen-Olesen vortex in the dual Ginzburg-Landau model and applying the universal 
topological invariants of the knot geometry being the length-to-diameter ratios of the fluxtubes. 
A number of meson fj states are identified as the knot-like gluonic excitations in the form of 
closed chromo-electric fluxtubes with the chromo-flux quantum n = 1 and the vibrationally- 
excited string mode m = 2. The most of the low-lying fj states are found to be described 
by the knotted string excitations of the types 2f and K = 3^ topologies, having the quantum 
number (n, m) = j = 0, 1, 2. 

Given that the glueballs can mix with the quark-antiquark states, as proposed in |41] . our 
calculation present a support that some of the (7^-states(mainly fj states and a few r/), listed 
in Table II and III, can have a dominate component of glues, namely, the energy of the valence 
quarks is negligible, and as a result, these fj states can be well described by the knotted fluxtubes 
formed by the collective gluons. The origin of the uncertainty in the knot-type identifications for 
a few states, such as /o(980) and /o(1370), /2(2220), r/(2225) in Table II and III, remains unclear 
yet, and the further studies are needed for the origin though the quark hybrid component of the 
states is probable origin. 

D. Jia is grateful to X. Liu for discussions. This work is supported in part by National 
Natural Science Foundation of China (No. 10965005), The Project-sponsored by SRF for ROCS, 
SEM, and by The Project of Key Laboratory of Atomic and Molecular Physics & Functional 
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